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We  have  studied  the  reactions  of  Ar\  02+,  and  CyH7+  with  toluene  (CH3-C6H5)  and  ethylbenzene 
(C2H5-C6H5)  as  a  function  of  temperature.  Ar+  reacts  with  C7H8  at  the  collision  rate  constant  at  both 
300  and  500  K  to  produce  mainly  C7H7+  cations  ( — 85%),  with  the  remainder  of  the  reactivity  divided 
into  numerous  small  channels.  The  fraction  of  C7H7+  produced  as  benzylium+  (6  membered  ring,  Bz+) 
as  opposed  to  tropylium+  (7  membered  ring,  Tr+)  is  determined  through  ion-molecule  reactivity  to  be 
primarily  Bz+.  The  reaction  of  Ar+  with  C8H10  proceeds  at  slightly  lower  than  the  collision  rate  constant 
at  300  K,  but  increases  to  the  collision  rate  constant  at  500  K.  Over  80%  of  the  products  were  C7H7+,  the 
large  majority  of  which  also  is  Bz+.  The  02+  reaction  with  CgHio  proceeds  at  the  collision  rate  constant 
and  produces  C8Hi0+,  CeH6+,  and  C7H7+,  with  the  latter  again  primarily  being  Bz+.  C7H7+  was  found  to 
react  with  both  CH3-C6H5  and  C2H5-C6H5  by  both  collisionally  stabilized  clustering  and  a  bimolecular 
reaction  forming  R-C6H4-CH2+.  At  300  K  the  rate  constants  are  about  half  of  the  collision  rate  constant 
and  drop  rapidly  with  temperature.  About  2/3  of  the  products  are  clusters  at  300  K,  1/3  at  400  K,  and 
negligible  amounts  at  500  and  600  K.  Evidence  for  thermal  dissociation  of  the  cluster  is  presented.  Calcu¬ 
lations  show  the  most  stable  structure  is  the  one  where  the  methylene  group  of  benzylium  is  coordinated 
at  the  para  position  of  the  alkylbenzene,  with  one  phenyl  group  approximately  45°  out  of  plane  relative 
to  the  other.  A  bond  strength  of  15-16kcalmol_1  is  calculated  for  both  neutrals.  Numerous  less  stable 
structures  were  found. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

It  has  long  been  known  that  C7H7+  exists  as  both  six  and  seven 
member  ring  structures  -  benzylium  (Bz+)  and  tropylium  (Tr+), 
respectively  [1-11].  These  rather  stable  ions  are  formed  in  many 
ion-molecule  reactions  of  alkylbenzenes,  and  the  history  of  their 
study  has  been  detailed  elsewhere  [1-13].  Bz+  reacts  with  alkyl¬ 
benzenes  by  reaction  (1): 

CH2C6H5+  +  R-C6H5^  CH2C6H4-R+  +  C6H6  (la) 

CH2C6H5+  +  R-C6H5  +  ]VI  -*  CH2C6H5(R-C6H5)+  +  M  (lb) 

while  Tr+  is  unreactive,  providing  a  convenient  way  to  determine 
how  much  of  each  isomer  is  present. 

We  have  recently  observed  C7H7+  while  developing  a  new 
flowing  afterglow  Langmuir  probe  technique  to  study  dissociative 
recombination  (i.e.  AB+  +  e_  ^  A  +  B)  [14].  This  is  a  variation  of  the 
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variable  electron  and  neutral  density  attachment  mass  spectrome¬ 
try  (VENDAMS)  approach  to  studying  a  wide  variety  of  previously 
inaccessible  areas  of  plasma  kinetics  [15].  In  this  technique,  Ar+  is 
allowed  to  react  with  a  neutral  (e.g.  toluene)  to  form  one  or  more 
polyatomic  ions.  The  ratios  of  Ar+  to  product  cations  are  monitored 
as  a  function  of  the  initial  plasma  density.  Because  Ar+  does  not 
recombine  with  e-,  but  the  polyatomic  ions  do,  at  higher  plasma 
densities  the  Ar+  to  product  ion  ratio  increases;  analysis  of  the 
shape  of  the  curve  describing  this  behavior  yields  the  rate  constants 
for  dissociative  recombination  of  the  various  cations.  In  our  initial 
paper  on  dissociative  recombination  [14  ,  we  chose  toluene  (C7H8) 
as  one  of  the  neutrals,  which  produces  C7H7+  as  the  major  product. 
However,  in  those  experiments  we  were  unable  to  determine  the 
ratio  of  Bz+/Tr+.  In  the  present  work  we  have  examined  the  reaction 
of  Ar+  with  toluene  as  a  function  of  temperature  in  order  to  deter¬ 
mine  the  Bz+/Tr+  product  ratios.  We  have  also  expanded  the  work 
to  include  ethylbenzene,  in  part  for  future  dissociative  recombina¬ 
tion  studies.  As  a  check,  we  have  studied  02+  with  ethylbenzene 
to  compare  to  previous  experiments  [16].  Further,  we  have  studied 
the  Bz+  reaction  with  both  toluene  and  ethylbenzene,  as  our  mod¬ 
eling  relies  on  this  secondary  chemistry.  Since  Bz+  clusters  to  both 
toluene  and  ethylbenzene  at  the  pressures  used  in  our  experiments, 
calculations  of  the  cluster  bond  strengths  for  reaction  (lb)  have 
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also  been  made.  Considerable  insight  into  the  reaction  mechanism 
is  obtained  from  the  measurements,  calculations,  and  comparisons 
to  previous  data. 


2.  Experimental  and  theoretical  methods 

The  experiments  were  performed  in  the  AFRL  variable  temper¬ 
ature  selected  ion  flow  tube  (SIFT).  The  instrument  has  been  well 
described  previously  [17],  and  only  details  important  to  the  present 
measurements  are  given.  Ions  are  made  in  a  high  pressure  ion 
source  through  electron  impact.  Pure  Ar  and  02  were  used  to  make 
Ar+  and  02+,  respectively.  A  mixture  of  Ar  and  C7H8  was  used  to 
make  C7H7+  in  the  source,  which  was  determined  to  be  20-30%  Tr+ 
with  the  remainder  Bz+.  The  ions  are  extracted  and  injected  into 
a  quadrupole  mass  filter  and  the  desired  mass  ion  was  selected. 
The  ions  enter  the  flow  tube  through  a  Venturi  inlet  and  are  car¬ 
ried  downstream  by  a  helium  buffer.  After  several  milliseconds  and 
approximately  104-105  collisions  with  helium,  C7H8  or  CgHio  was 
added  through  a  finger  inlet.  Reactions  occurred  over  the  remaining 
59  cm  of  the  flow  tube,  and  subsequently  the  flow  was  sampled 
through  a  small  orifice  in  a  rounded  nose  cone.  The  remainder  of 
the  flow  was  pumped  away  through  a  MKS  servo-controlled  gate 
valve,  which  maintained  a  controlled  pressure  within  the  flow  tube. 
Sampled  ions  were  then  injected  into  a  second  quadrupole  mass 
spectrometer  and  detected  by  a  discrete  dynode  particle  multiplier. 
The  entire  flow  tube  can  be  heated  by  resistance  heaters  in  several 
zones.  Here,  most  of  the  measurements  were  only  done  at  300  K  and 
500  K  because  the  temperature  dependence  is  weak.  The  exception 
is  the  work  on  C7H7+,  for  which  the  temperature  dependences  are 
strong. 

All  flows  were  measured  using  MKS  mass  flow  controllers.  The 
ion  velocity  was  calculated  from  the  helium  flow  velocity  and  a  pre¬ 
viously  measured  ion  to  helium  velocity  ratio  of  ~1.6  [18 1.  Coupled 
with  the  known  reaction  length  of  59  cm,  this  yields  the  reaction 
time.  Rate  constants  were  obtained  by  monitoring  the  primary 
decay  as  a  function  of  added  neutral  concentration  in  the  flow  tube. 
For  Ar+  and  02+,  single  exponential  decays  were  found  and  rate 
constants  obtained  in  the  normal  fashion.  For  C7H7+,  the  decays 
were  fit  to  an  exponential  plus  a  constant  functional  form.  The 
decaying  portion  was  due  to  the  Bz+  reactions  and  the  constant 
was  attributed  to  Tr+. 

Toluene  and  ethylbenzene  were  introduced  both  neat  and  as 
mixtures  in  helium.  The  total  count  rate  was  largely  independent 
of  the  extent  of  reaction,  indicating  minimal  mass  discrimination 
in  the  data,  and  no  correction  was  needed.  Branching  ratios  (with¬ 
out  distinguishing  between  C7H7+  isomers)  were  determined  from 
the  early  part  of  the  decay  in  the  normal  SIFT  manner  of  extrap¬ 
olating  product  branching  fractions  to  zero  flow.  Anderson  et  al. 
showed  that  as  long  as  the  amount  of  primary  ion  depletion  was 
kept  small,  the  technique  yields  nascent  branching  ratios  to  within 
a  few  percentage  points  [19].  Alternatively,  we  iteratively  solve  the 
set  of  coupled  differential  equations  describing  the  reaction  scheme 
throughout  the  known  reaction  time,  comparing  the  calculated 
product  ion  abundances  to  those  measured  through  a  weighted 
least  squares  fit  [15].  By  varying  all  relevant  rate  constants  through 
a  Monte  Carlo  procedure,  over  ranges  constrained  either  by  the 
literature  or  calculated  collision  rate  constants,  the  full  parame¬ 
ter  space  may  be  efficiently  sampled.  The  latter  analysis  was  done 
mainly  to  determine  the  Bz+/Tr+  ratio.  The  best  fit  value  of  the  rate 
constant  for  each  reaction  channel  is  determined  by  identifying 
the  minimum  in  the  least  squares  parameter  as  a  function  of  that 
rate  constant;  uncertainties  are  determined  by  values  at  the  largest 
least  squares  parameter  that  still  provide  an  acceptable  description 
of  the  data.  We  find  good  agreement  in  product  branching  ratios 
determined  using  the  two  methods. 


Errors  for  the  linear  decay  rate  constants  are  ±25%  for  absolute 
values  and  ±15%  for  relative  changes  as  a  function  of  temperature. 
The  error  in  branching  is  estimated  to  be  5-1 0  percentage  points  for 
the  large  channels  and  1  -2  percentage  points  for  the  small  chan¬ 
nels.  For  the  Bz+  reactions  errors  are  larger  and  are  discussed  in 
Section  3.4.  We  estimate  that  the  Bz+  and  Tr+  branching  are  accurate 
to  10  percentage  points. 

Calculations  of  the  Bz+-toluene  and  Bz+-ethylbenzene  cluster 
bond  strengths  were  performed  within  the  Gaussian  09W  package 
at  the  B3LYP/6-31 1  +G(d)  level  [20-22].  Several  starting  geome¬ 
tries  were  explored,  including  geometries  found  in  the  literature 
for  similar  cluster  ions  [23-27].  Each  optimized  structure  was  con¬ 
firmed  as  a  true  minimum  via  a  frequency  analysis.  All  clusters 
possess  singlet  ground  state  structures,  and  energies  we  provide 
include  zero-point  corrections. 


3.  Results  and  discussion 

3.1.  Reactions  of  02+  with  C8HW 

As  a  test  of  the  system,  we  have  repeated  our  previous  work  on 
the  reaction  of  02+  with  C8H10.  The  results  are  compared  in  Table  1 . 
The  rate  constants  are  near  the  collision  rate  constant  (within 
~1 0%)  and  in  good  agreement  with  several  previous  measurements 
[16,28,29].  The  branching  measurements  are  generally  in  good 
agreement  with  previous  measurements  [6,16,28,29]  with  about 
75%  of  the  reactions  producing  C7H7+,  the  bulk  of  the  remainder 
being  non-dissociative  charge  transfer,  and  finally  2-3%  of  C6H6+. 
No  clear  trend  with  temperature  is  apparent;  the  small  changes 
are  within  our  estimated  error.  Less  non-dissociative  charge  trans¬ 
fer  was  found  at  high  temperature,  as  one  would  expect,  but  the 
observed  change  is  within  our  error.  The  division  of  C7H7+  into  Bz+ 
and  Tr+  is  in  excellent  agreement  with  the  ICR  values  in  Fridgen 
et  al.  [6]  and  reasonable  agreement  with  the  SIFT  values  in  that 
paper.  Poorer  agreement  is  found  with  the  ratio  reported  by  Arnold 
et  al.  [28].  The  present  measurements  employed  different  methods 
of  extracting  the  Bz+/Tr+  ratio  than  that  used  in  the  previous  SIFT 
experiments.  The  present  measurement  and  the  ICR  values  of  Frid¬ 
gen  et  al.  have  a  lower  contribution  from  Tr+  than  the  previous  SIFT 
measurements.  The  agreement  between  the  current  study  and  Frid¬ 
gen  et  al.  [6],  as  well  as  the  detailed  statistical  analysis  performed 
here,  lead  us  to  favor  the  present  results. 


3.2.  Reactions  of  Ar+  with  C7H8 

Our  results  for  the  reaction  of  Ar+  with  C7H8  are  listed  in  Table  2. 
The  overall  rate  constant  is  within  6%  of  the  collision  rate  constant 
[30]  at  both  300  and  500  K,  well  within  our  uncertainty.  The  vast 
majority  of  products  are  the  two  C7H7+  cations  (~85%).  All  product 
ions  can  be  described  as  charge  transfer,  most  of  it  dissociative;  only 
about  2%  of  the  non-dissociative  channel  is  left  after  correcting  the 
92  amu  peak  for  13C  containing  C7H7+.  Products  with  4,  5,  6,  and  7 
carbons  were  found.  The  largest  non-C7H7+  cation  is  C5H6+  at  3-4%. 
For  reactions  with  known  thermochemistry,  the  exothermicity  is 
64kcalmol_1  and  higher,  so  it  is  not  surprising  that  dissociative 
charge  transfer  is  the  main  process.  About  2-3  times  more  Bz+  is 
formed  than  Tr+,  with  the  former  accounting  for  about  60%  of  the 
products  at  both  temperatures  and  the  latter  25%.  While  there  is 
some  scatter  between  the  temperatures,  it  is  hard  to  discern  trends 
because  many  of  the  products  are  so  small  and  their  dependences 
are  weak.  Considering  the  exothermicity  of  the  reaction,  it  is  not 
surprising  that  relatively  small  amounts  of  added  thermal  energy 
have  only  a  small  effect. 
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Table  1 

Product  distributions  and  exothermicities  for  the  reaction  of  02+  with  C8Hi0  (ethylbenzene).  Product  distributions  are  corrected  for  isotopes,  but  not  corrected  for  impurity 
ions  from  H20+  or  air  related  ions.  The  measured  total  rate  constant  and  calculated  collision  rate  constant  [30]  are  shown  at  the  bottom  of  the  table.  The  error  in  branching  is 
estimated  to  be  5-10  percentage  points  for  the  large  channels  and  1-2  percentage  points  for  the  small  channels;  measured  rate  constants  have  error  bars  of  ±25%.  Literature 
values  are  shown  in  brackets. 


Products 

Exothermicity  (kcal  mol-1 ) 

300  K 

%  branching  [lit.  values] 

500  K 

%  branching  [lit.  values] 

c6h6+ 

39.7 

3  [3]a,b,c 

2  [2.5]b 

C7H7+  (total) 

71  [67]a,b  [61  ]c 

76  [75.5]b 

Bz+ 

59.0 

63  [44]a  [50]  [64]d 

71 

Tr+ 

48.0 

8  [23]a  [11]  [7]d 

5 

c8h10+ 

76.3 

25  [30]a,b  [36]c 

[20]b 

Total  rate  constant  (cm3  s-1 )  [lit.  values] 

1.7  x  10“9[1.7  x  1 0— 9 ]a[2.2  x  10-9]a-b 

1.6  x  10“9  [2.2  x  10_9]b 

Collision  rate  constant  (cm3  s_1 ) 

1.9  xl0“9 

1.9  x  10“9 

a  Arnold  et  al.  [28]. 

b  Williams  et  al.  [16]  (average  of  two  flow  tube  values,  when  appropriate). 
c  Fridgen  et  al.  [6]  SIFT  measurement  at  1  Torr. 

d  Fridgen  et  al.  [6]  ICR  measurement,  assumes  total  C7H7+  is  the  current  value  of  71%. 


Table  2 

Product  distributions  and  exothermicities  for  the  reaction  of  Ar+  with  C7H8 
(toluene).  Product  distributions  are  corrected  for  isotopes,  but  not  corrected  for 
impurity  ions  from  H20+  or  air  related  ions.  The  measured  total  rate  constant  and 
calculated  collision  rate  constant  [30]  are  shown  at  the  bottom  of  the  table.  The  error 
in  branching  is  estimated  to  be  5-10  percentage  points  for  the  large  channels  and 
1-2  percentage  points  for  the  small  channels;  measured  rate  constants  have  error 
bars  of  ±25%. 


Products 

Exothermicity 
(kcal  mol-1 ) 

300  K 

%  branching 

500  K 

%  branching 

c4h4+ 

0.6 

1.5 

c5h4+ 

0.2 

0.4 

c5h5+ 

64.1 

1.5 

2.5 

c5h6+ 

89.9 

3.4 

4.1 

c6h5+ 

2.1 

3.2 

c7h5+ 

0.4 

1.1 

c7h6+ 

100.5 

2.0 

2.0 

C7H7+  (total) 

87.6 

83.0 

Bz+ 

108.6 

64 

56 

Tr+ 

97.5 

23 

27 

c7h8+ 

160.0 

2.0 

2.2 

Total  rate  constant 

1.5  xlO-9 

1.5  xlO”9 

(cm3  s_1) 

Collision  rate 

1.6  x  10“9 

1.6  x  10“9 

constant  (cm3  s_1) 

3.3.  Reactions  ofAr+  with  C8HW 

Table  3  shows  the  results  for  Ar+  with  C8Hi0.  The  rate  constants 
are  near  the  collision  rate  constant  at  500 1<  (95%),  but  not  quite  at 
the  collision  rate  constant  at  300  K  (72%).  This  difference  is  outside 
our  relative  error  and  the  unusual  positive  temperature  depend¬ 
ence  for  a  highly  exothermic  ion-molecule  reaction  was  confirmed 


through  several  repeated  measurements.  Possibly  the  large  energy 
mismatch  between  the  Ar  and  C8Hi0  ionization  energies  is  respon¬ 
sible,  resulting  in  curve  crossings  that  are  not  completely  favorable 
to  charge  transfer.  The  reaction  proceeds  entirely  through  dissocia¬ 
tive  charge  transfer.  This  is  in  contrast  to  the  02+  reaction  where  a 
quarter  of  the  reaction  was  non-dissociative.  Again  the  large  recom¬ 
bination  energy  of  Ar+  is  almost  certainly  the  reason.  Just  over  80% 
of  the  products  are  C7H7+  at  both  temperatures.  The  remainder  of 
the  products  are  split  nearly  equally  between  C6H6+  and  C8H9+.  At 
300  K,  we  find  the  Bz+  to  Tr+  ratio  to  be  about  8,  in  good  agreement 
with  Fridgen  et  al.  [6]  There  was  not  enough  curvature  in  the  data  at 
500  K  to  obtain  the  Bz+/Tr+  ratio.  However,  based  on  the  02+  results, 
we  do  not  expect  much  change  from  the  300  K  result. 

3.4.  Reactions  ofBz+  with  C7H8  and  C8HW 

The  final  reactions  studied  were  Bz+  with  both  C7H8  and  C8Hi0. 
Both  reactions  involve  the  scheme  detailed  in  Section  1,  namely 
R  and  H  transfer  (la)  and  association  (lb).  As  shown  later,  the 
Bz+(R-C6H5)  cluster  bond  strengths  are  rather  small,  indicating 
that  the  clusters  are  barely  stable  with  respect  to  thermal  disso¬ 
ciation  at  room  temperature,  and  become  increasingly  unstable  as 
the  temperature  is  raised.  Therefore,  the  complete  reaction  mech¬ 
anism  also  involves  the  reverse  of  (lb),  which  is  hidden  in  the  data. 
Additionally,  in  an  attempt  to  measure  pressure  dependences,  it 
was  found  that  the  overall  rate  constants  decreased  slightly  with 
pressure.  This  is  not  physical  and  indicates  that  the  fragile  species 
were  partially  dissociating  due  to  the  electric  fields  used  in  samp¬ 
ling.  For  these  reasons,  our  uncertainties  are  larger  than  normal  - 
we  estimate  the  data  are  only  accurate  to  a  factor  of  two.  However, 
the  trends  are  so  large  the  data  are  still  informative. 


Table  3 

Product  distributions  and  exothermicities  for  the  reaction  of  Ar+  with  C8Hi0  (ethylbenzene).  Product  distributions  are  corrected  for  isotopes,  but  not  corrected  for  impurity 
ions  from  H20+  or  air  related  ions.  The  measured  total  rate  constant  and  calculated  collision  rate  constant  [30]  are  shown  at  the  bottom  of  the  table.  The  error  in  branching  is 
estimated  to  be  5-10  percentage  points  for  the  large  channels  and  1-2  percentage  points  for  the  small  channels;  measured  rate  constants  have  error  bars  of  ±25%.  Literature 
values  are  shown  in  brackets. 


Products 

Exothermicity  (kcal  mol-1 ) 

300  K 

%  branching 

500  K 

%  branching 

c6h6+ 

124.8 

8 

8 

C7H7+  (total) 

81 

84 

Bz+ 

167.2 

72  [68]a 

b 

Tr+ 

156.1 

9  [13]a 

b 

c8h9+ 

119.9 

11 

7 

c8h10+ 

161.4 

c 

c 

Total  rate  constant  (cm3  s-1 )  [lit.  values] 

1.3  x  10_9[1.9  x  10_9]a 

1.7  x  10“9 

Collision  rate  constant  (cm3  s_1 ) 

1.8  x  10-9 

1.8  x  10-9 

a  Fridgen  et  al.  [6]  ICR  measurement,  assumes  total  C7H7+  is  the  current  value. 
b  Not  enough  curvature  for  accurate  determination. 
c  Too  small  to  measure  after  correcting  for  13CC7Hg+. 
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300  K 


Fig.  1.  Plot  of  C7H7+  signal  as  a  function  of  toluene  concentration  at  300  K.  The 
fraction  of  C7H7+  that  declines  is  due  to  benzylium  and  the  flat  portion  is  due  to 
tropylium. 


Table  4 

Rate  constants  and  product  distributions  for  the  reaction  of  Bz+  with  C7H8  and  C8Hi0 
at  a  number  density  of  3  x  1016  atoms  cm-3  helium.  The  cluster  product  is  shown  in 
the  table,  the  reactive  fraction  is  the  remainder.  Branching  and  rate  constants  are 
accurate  to  a  factor  of  two  (see  Section  3.4). 


Temp  (I<) 

c7h8 

c8h10 

k  (cm3  s_1) 

Cluster  (%) 

k  (cm3  s_1 ) 

Cluster  (%) 

300 

7  x  10-10 

60 

6  x  10-10 

68 

400 

4  x  10-10 

31 

4  x  10-10 

30 

500 

1  xlO-10 

1 

1  x  10-10 

1 

600 

3  x  10-11 

0 

5  x  10-11 

1 

Fig.  1  shows  a  decay  of  C7H7+  as  a  function  of  C7H8  concentra¬ 
tion  at  300  K.  The  C7H7+  initially  declines  and  then  levels  out.  The 
decline  is  due  to  reaction  of  Bz+  and  the  flat  portion  is  Tr+.  Typically 
about  1/3  of  the  ions  were  Tr+.  The  rate  constants  are  derived  by 
fitting  to  an  exponential  plus  a  constant,  and  the  line  in  Fig.  1  is 
the  least  squares  fit.  Because  only  a  single  species  reacts,  product 
branching  can  be  deduced  in  the  normal  manner.  Table  4  shows 
the  rate  constants  and  branching  percentages  for  Bz+  reacting  with 
C7FI8  and  C8H10  at  a  number  density  of  3  x  1016  atoms  cm-3  helium. 
Overall  rate  constants  decline  sharply  from  300  to  600  K,  by  fac¬ 
tors  of  23  and  12  for  toluene  and  ethylbenzene,  respectively.  At 
300  K,  60%  and  68%  of  the  reactivity  with  toluene  and  ethylben¬ 
zene,  respectively,  produces  the  cluster  and  the  rate  constants  are 
a  large  fraction  of  the  collision  rate  constant  (~l/3).  By  400  K,  the 
rate  constants  drop  by  about  50%,  as  does  the  amount  of  clustering 
product.  There  is  essentially  no  cluster  left  at  500  and  600 1<  (<1%) 
and  much  greater  reductions  in  the  rate  constants  are  found.  The 
kinetics  for  the  two  reactions  are  similar,  which  is  not  surprising 
given  the  similarity  between  the  neutral  species.  Thermal  dissoci¬ 
ation  would  normally  show  as  curved  decay  plots.  Since  the  data 
already  have  curvature  due  to  the  two  isomers,  it  is  not  possible  to 
quantify  a  second  source  of  curvature  and  no  information  on  the 
thermal  dissociation  rates  was  derived. 


There  have  been  several  previous  studies  of  the  kinetics  of  Bz+ 
reacting  with  C7H8,  and  one  of  Bz+  reacting  with  C8Hi0.  Dunbar 
[4]  found  only  reaction  la  at  1.8  x  10-8  Torr,  and  determined  a  rate 
constant  of  4  x  10-10  cm3  s-1  for  C7FI8.  Ausloos  et  al.  [10]  studied 
both  the  C7IT8  and  CsF^o  reactions  at  low  pressure.  They  found  val¬ 
ues  of  1 .7  x  1 0-10  cm3  s-1  and  5.7  x  1 0-10  cm3  s-1 ,  respectively.  In 
agreement  with  Dunbar,  they  observed  only  CFI2  transfer  for  the 
C7IT8  reaction.  Giardini-Guidoni  and  Zocchi  [11]  found  a  very  small 
rate  constant  in  a  radiolysis  experiment  (2.1  x  10-11  cm3  s-1 );  this 
value  is  clearly  an  outlier.  For  the  CsF^o  reaction,  they  found 
that  ~3/4  of  the  reactions  proceed  by  CH2  transfer  and  1/4  by 
a  hydride  transfer.  For  the  C7FI8  reaction,  Sharma  and  Kebarle 
[7]  found  only  clustering  at  300 K  and  4Torr  CFI4,  and  a  faster 
rate  constant  of  2  x  10-9  cm3  s-1 .  At  436  K,  they  observed  that  the 
rate  constant  decreased  substantially  and  that  thermal  dissocia¬ 
tion  occurred  without  any  quantification  given.  Combining  all  the 
C7IT8  data  shows  that  the  reaction  is  slower  at  lower  pressures  and 
produces  only  the  CFI2  transfer  product.  In  contrast,  the  rate  con¬ 
stant  increases  to  the  collision  rate  constant  at  4 Torr  (CH4)  and 
produces  only  the  cluster.  The  present  result  at  0.9  Torr  helium  is 
intermediate  in  both  rate  constant  and  branching.  We  did  not  see 
the  hydride  transfer  channel  for  the  CsFIto  reaction. 

In  order  to  understand  the  data,  we  have  calculated  the  bond 
strengths  of  both  Bz+-C7H8  and  Bz+-C8FIio.  Numerous  stable  struc¬ 
tures  were  found  and  Fig.  2  shows  the  structures  and  calculated 
stabilities.  The  most  stable  structure,  shown  at  the  bottom  of 
the  figure,  is  similar  for  both  clusters.  The  methylene  group  of 
the  benzylium  structure  is  coordinated  at  the  para  position  of 
the  alkylbenzene,  with  one  phenyl  group  approximately  45°  out 
of  plane  relative  to  the  other.  We  find  the  bond  strengths  for 
these  structures  are  15.2  and  15.8  kcal  mol-1  for  Bz+-toluene  and 
Bz+-ethylbenzene,  respectively.  We  refer  to  those  as  “para”  in 
the  figure.  For  C7H8,  we  found  an  “ortho”  structure  where  the 
methylene  of  benzylium  is  coordinated  at  the  ortho  position  of  the 
alkylbenzene.  This  is  only  a  couple  of  kcal  mol-1  less  stable  than  the 
“para”  structure.  No  similar  structure  was  found  for  C8Hi0,  although 
it  was  searched  for  extensively. 

The  least  stable  structures  found  are  the  “stacked-same”  struc¬ 
ture  (~5.5  kcal  mol-1 )  for  both  neutrals  and  the  “stacked-opposite” 
structure  for  C7H8.  These  structures  have  the  phenyl  rings  stacked, 
and  the  functional  groups  are  on  similar  carbons  (“same”),  or  oppo¬ 
site  carbons  (“opposite”),  in  their  respective  phenyl  rings.  The 
stacked  opposite  structure  for  C8Hi0  is  over  twice  as  stable  as 
that  for  C7H8.  A  similarly  stable  “stacked-opposite  offset”  structure 
was  found  for  C8Hi0.  “T-shape”  structures,  where  the  benzylium  is 
slightly  off-vertical  relative  to  the  horizontal  alkylbenzene  phenyl 
ring,  were  found  to  be  almost  as  weakly  bonded  as  the  stacked- 
same  structures  (~5.8  kcal  mol-1 ). 

In  our  flow  tube  only  the  most  stable  “para”  structures  are  sta¬ 
ble  enough  to  survive  the  thousands  of  helium  collisions  [31,32]. 
To  emphasize  this  point,  we  have  calculated  the  energy  distribu¬ 
tions  of  the  most  stable  Bz+-C7H8  cluster  at  various  temperatures. 
The  results  are  shown  in  Fig.  3.  The  dashed  vertical  line  represents 
the  bond  dissociation  energy.  At  300  K,  most  Bz+-C7FI8  molecules 
have  energies  below  the  dissociation  energy  (99%).  Given  the  multi¬ 
collision  effects  in  thermal  dissociation,  even  at  this  temperature 
some  thermal  dissociation  may  be  occurring.  At  400  K,  30%  of  the 
molecules  have  enough  energy  to  dissociate.  While  the  reactive 
channel  is  roughly  constant  from  300  to  400 1<  (partial  rate  constant 
~2.8  x  10-1°  cm3  s-1),  the  partial  rate  constant  for  the  clustering 
channel  decreases  by  about  a  factor  of  4.  Some  of  that  decrease  is 
probably  related  to  thermal  decomposition  affecting  the  decline  in 
the  Bz+.  At  500  K,  there  is  essentially  no  cluster  left;  the  partial  rate 
constant  drops  by  two  orders  of  magnitude  and  79%  of  the  cluster 
has  enough  energy  to  dissociate.  The  reactive  channel  also  drops 
by  a  factor  of  three.  The  trends  continue  at  600 1<:  no  association 
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Fig.  2.  Optimized  geometries  of  the  isomers  for  the  benzylium-toluene  (C7H8-C7H7+,  left)  and  benzylium-ethylbenzene  (C8Hio-C7H7+  right)  clusters.  Carbon  atoms  and 
hydrogen  atoms  are  gray  and  white,  respectively.  Numbers  are  the  bond  strengths  in  kcal  mol-1. 


product  was  found  and  the  reactive  channel  decreased  by  another 
factor  of  three.  The  steepness  of  the  temperature  dependence  of 
the  clustering  channel  is  much  larger  than  can  be  expected  without 
consideration  of  thermal  dissociation. 

The  400-600  temperature  dependence  of  the  reactive  channel 
is  also  large,  suggesting  a  more  complicated  mechanism  involving 
thermal  activation.  An  example  of  this  has  been  seen  by  our  group 
previously  33].  At  300 K,  C1_(H20)  reacts  with  CH3Br  to  produce 
almost  exclusively  Br-.  That  is  an  endothermic  channel.  At  lower 
temperatures,  Cl~(CH3Br)  was  observed,  showing  that  the  reaction 
proceeded  by  a  very  slightly  endothermic  ligand  switching  to  form 
Cl-(CH3Br),  which  then  underwent  thermal  activation  to  form  Br¬ 
in  competition  with  Cl-  formation. 

We  speculate  that  a  similar  mechanism  occurs  here.  The  com¬ 
plete  mechanism  would  then  include  the  following  steps.  First,  a 
complex  is  made. 

+CH2-C6H5  +  R-C6H5  -*  (CH2-C6H5-R-C6H5)+*  (3) 


The  fact  that  4Torr  of  CH4  is  enough  to  completely  stabilize  the 
complex  at  300  K  and  0.9Torr  Helium  only  partially  stabilizes  it 
indicates  that  the  complex  lifetime  is  on  the  order  of  the  helium 
collision  time  (i.e.  about  10-7  s).  The  reactive  channel  must  also  be 
this  slow  or  it  would  occur  before  the  complex  can  be  stabilized. 
The  excited  complex  can  then  either  dissociate  (-3),  go  to  products 
via  transfer  of  the  R-group  (4),  or  be  stabilized  (5). 

(CH2-C6H5-R-C6H5)+*  ->  +CH2-C6H5  +  R-C6H5  (-3) 

(CH2-C6H5-R-C6H5)+*  ->  +CH2-C6H4-R  +  C6H6  (4) 

(CH2-C6H5-R-C6H5)+*  +  He  ->  (CH2-C6H5-R-C6H5)+ +  He  (5) 

Additionally,  the  stable  complex  can  be  activated,  particularly 
at  high  temperatures 

(CH2-C6H5-R-C6H5)+  +  He  -*  (CH2-C6H5-R-C6H5)+*  +  He  (-5) 
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Internal  energy  (eV) 


Fig.  3.  Calculated  thermal  energy  distributions  of  the  Bz+-C7H8  “para”  cluster  at 
the  indicated  temperatures.  The  bond  dissociation  energy  (BDE,  15.2  kcal  mol-1 )  is 
indicated  by  the  vertical  line,  and  the  fraction  of  each  distribution  above  the  BDE  is 
indicated  in  parentheses. 

At  that  point,  reactions  (-3),  (4),  and  (5)  occur.  The  thermally 
activated  reactive  channels,  reaction  (-5)  followed  by  reaction  (4), 
are  calculated  to  be  7.8  and  7.6  kcal  mol-1  endothermic  for  the  most 
stable  isomers.  These  are  considerably  lower  in  energy  than  the 
simple  cluster  dissociations  represented  by  reactions  (-5)  followed 
by  (-3 ),  which  are  about  1 5  kcal  mol-1  endothermic.  A  considerably 
high  transition  state  barrier  along  the  path  of  (-5),  (4)  could  alter  the 
conclusion,  but  locating  that  stationary  point  is  beyond  the  scope 
of  this  work. 

4.  Conclusions 

We  have  explored  several  aspects  of  the  reactivity  involving 
toluene,  ethylbenzene,  and  C7H7+.  Reactions  of  Ar+  with  toluene 
and  ethylbenzene  have  been  probed  to  determine  the  ratio  of  ben- 
zylium  to  tropylium.  The  results  show  that  mainly  Bz+  is  formed, 
in  addition  numerous  minor  products  were  observed.  Very  lit¬ 
tle  dependence  on  temperature  was  found,  probably  because  the 
charge  transfer  energy  is  so  large.  Similarly,  we  explored  the  02+ 
reaction  with  ethylbenzene  and  helped  solve  a  previous  discrep¬ 
ancy  in  the  amount  of  Tr+  formed.  Finally,  we  explored  temperature 
dependences  of  the  reactions  of  Bz+  with  toluene  and  ethylbenzene. 
Combining  our  results  at  0.9Torr  helium  with  previous  studies  at 
both  higher  and  lower  pressure  shows  that  the  reaction  is  complex, 
involving  clustering,  R-group  transfer,  and  thermal  activation. 
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